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SUMMARY 


A relatively simple equation is presented for estimating the induced drag 
ideal efficiency factor e for arbitrary cross-sectional wing forms. This 
equation is based on eight basic but varied wing configurations which have 
exact solutions. The e function which relates the basic wings is developed 
statistically and is a continuous function of configuration geometry. The 
basic wing configurations include boxwings shaped as a rectangle, ellipse, and 
diamond; the V-wing; end-plate wing; 90 degree cruciform; circle dumbbell; and 
biplane. Example applications of the e equations are made to many wing forms 
such as wings with struts which form partial span rectangle, diamond, triangle, 
ellipse, and semi -ell ipse wings; the ellipse and rectangle dumbbell wings; bow- 
tie, cruciform, winglet, and fan wings; and multi -wings. Derivations are 
presented in the appendices of exact closed form solutions found of e for the 
V-wing and 90 degree cruciform wing and for an asymptotic solution for multi- 
wings . 


INTRODUCTION 

With the objective of better economics through improved aerodynamic and 
structural shapes, many unconventional aircraft configurations are being 
investigated. Particularly desirable are wing forms for which both drag and 
structural stresses decrease. A measure of aerodynamics efficiency is the 
induced drag ideal efficiency factor denoted by e. 

Historically, exact solutions for induced drag have been found for several 
nonplanar wing configurations. These include the rectangle boxwing, end-plate 
wing, and biplane reported in reference 1; the ellipse boxwing and circle 
dumbbell wing from the methods of reference 2; and the diamond boxwing and 
V-wing in reference 3. These solutions form the basis for the present study. 

The objective of this study is to develop a generalized equation for esti- 
mating the ideal e for arbitrary yz-plane wing forms. For a given wing form 
the ideal e is the maximum that e can become, and at ideal e all the wing form 
surfaces have ideal loading distributions such that the total induced drag is 
minimized. A generalized equation for estimating e provides the designer with 



a tool for evaluating the aerodynamic efficiency for the configuration under 
consideration. 
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SYMBOLS 

lateral span along y-axis of wing configuration 
induced drag, equation (1) 

complete elliptic integral of the second-kind, with modulus k 
induced drag ideal efficiency factor, equation (1) 
e for end-plate wing, equation (16) 

e for V-wing, equation (15) 

function of parameter P, equation (20) 
function of n-|gH, equation (18) 

function of e , e , n equation (19) 

V G cl T 

dimensionless maximum height, h„/b 

m 

average dimensionless height, equation (25) 
maximum height of wing configuration 

complete elliptic integral of the first kind, with modulus k 
modulus for K and E, k = sina where a is the modular angle 

total lift of wing system, equation (1) 

total number of wings 

parameter, equation (21) 

perimeter of wing system, equation (29) 

perimeter of ellipse, equation (22) 

1 2 

free stream dynamic pressure, ^ pV 
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lateral coordinate, zero at midspan, positive to right 
gamma function, equation (14) 
dihedral angle 

dimensionless lateral coordinate, 2y/b 

averaged n of outboard face of wing, equations (28) and (76) 
span station of difference, equation (27) 
for fin wings, equations (30), (31), (32) 

a 

span stations, equation (27) 

span station of multi-wings, equation (76) 

air density 

perimeter parameter, equations (29) and (77) 

perimeter parameter of ellipse, equations (22) through (24) 

perimeter parameter of multi-wing, equation (78) 


IDEAL EFFICIENCY FOR BASIC CONFIGURATIONS 


The present analysis is based on a group of varied wing configurations 
which have exact solutions for predicting e. These include the wings described 
as rectangle boxwing, ellipse boxwing, diamond boxwing, V-wing, end-plate wing, 
90 degree cruciform wing, circle dumbbell wing, and biplane. The induced drag 
ideal efficiency factor listed throughout the paper relates the total induced 
drag to the total lift of the wing system by the equation: 



( 1 ) 


where b is the maximum span in the lateral plane, q is the dynamic pressure, and 
L is the sum of all the lifts, including wing, winglets, end-plates, struts, 
or tail. For ideal e, the lifts and loadings will have ideal distributions. 
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Rectangle and Ellipse Boxwings 

From reference 1, p. 205, for the rectangle boxwing, the eguations for 
the dimensionless maximum height for H and for e can be developed. 



E - (1 - k^)K 
E' - k^K' 


e 



( 2 ) 

(3) 

-b/2 0 b/2 



where K and E are complete elliptic integrals of the first and second kinds, 
respectively, with modulus k, while K' and E' are the complementary integrals. 
That is, K' is K with modulus A - k'2, and E' is E with modulus /I - k2. The 
modulus k is related to the modular angle a by k = sina. Eguations (2) and 
(3) do not permit a direct solution of e for a specified value of H. A value 
of k is specified and e from equation (3) is determined, then the corresponding 
H from equation (2) is determined using the same k value. Geometric character- 
istics of the rectangle include the nondimensional perimeter given by 


a 



2b + 2h 

m 

2b 


1 + H 


(4) 


The rectangle averaged cross-sectional area per wing span squared is 



The e for the ellipse boxwing can be 
determined from the work given in reference 
2. The results are 

e = 1 + H 


(5) 



( 6 ) 


The perimeter parameter of the ellipse is an elliptic integral of the second 
kind given by 


a 


e 


■ E, k = A - H2, Mil' 
HE, k = H > 1 


(7) 
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The ellipse averaged cross-sectional area per wing span squared is 


H 


a 


TTh b 
m 




(8) 


Relating Rectangle and Ellipse Boxwings 

Examination of equations (2) and (3) shows that at H = 0, e = 1, and at 
H 00, e 4H/ir. Then for the rectangle boxwing, e can be expressed as 

e = (1 + 7 H)(l + f) (9) 

where f is a function of H to be determined. Equation (9) can apply to both 
rectangle and ellipse boxwings by letting H = and f be a function of peri- 
meter difference. Then 


e = (1 + 7 H )(1 + f) 

TT a 

f = 1.4512 (1 - /P)^-°795 

P = a - Ce 



where is given in equation (7), and a for the rectangle is given in equation 
(4). The f function is statistically developed from the e and H data of equa- 
tions (2) and (3). Equation (10) applies to either the rectangle or ellipse 
boxwing and gives an estimate of e for any shape between these two. Example 
values of e are listed in table 1. 
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Then the dimensionless added mass is 
M 


1 

7T 




pira 


[r(^i) r(l -1)]‘ 


( 12 ) 


[r(^ + ^ tan * H) r(l - ^ tan * H)J^ it(1 - ^ tan"^ H) 


1 -1 .,vT2 


2 . 


where a is derived in reference 3, and H = h^^/b. In reference 3, D^- and L 
related by 


D,. 




^ 2MV^ 4a^ / M ^ 

irqb ( p) 

b pTia 

where from the sketch, 4a^/b^ = 1 + then 


e = 0 + H^) (M ) 


n(l + H^) (1 - - tan"' H) 

TT 


pira^ |^r(^ + ^ tan~^ H) r(l - ^ tan"^ H)j‘^ 


(13) 


(14) 


where r is the gamma function. For the diamond boxwing, values of e for 
various values of H are presented in table 1. Tables of the gamma function 
are given in reference 4. 


TABLE 1. - IDEAL e VALUES FOR RECTANGLE, ELLIPSE, 
AND DIAMOND BOXWINGS 


H 

rectangle 
eq. (10) 

el li pse 
eq. (6) 

di amond 
eq. (14) 

H 

rectangle 
eq. (10) 

ellipse 
eq. (6) 

di amond 
eq. (14) 

0 

1 

1 

— ^ 

1 

.8 

2.47997 

1.8 

1.29518 

.05 

1.15178 

1.05 

1.01270 

1.0 

2.78642 

2.0 

1.39320 

.10 

1 .26814 

1.10 

1.02622 

1.2 

3.08606 

2.2 

1 .49715 

.15 

1.37327 

1.15 

1.04056 

1 .4 

3.38041 

2.4 

1.60562 

.20 

1.47189 

1.20 

1.05573 

1.6 

3.67066 

2.6 

1 .71751 

.25 

1 .56610 

1.25 

1.07171 

1.8 

3.96086 

2.8 

1.83206 

.3 

1.65709 

1.3 

1 .08850 

2 

4.24189 

3 

1 .94866 

.4 

1.83204 

1.4 

1 .12442 

3 

5.63505 

4 

2.55051 

.5 

2.00026 

1.5 

1 .16330 

4 

6.99909 

5 

3.16840 

.6 

2.16362 

1.6 

1.20490 

oo 

-^4H/tt 

-41 

^2H/tt 


are 
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V-Wing 


An exact closed function solution was found for 
the e of a V-wing. The derivation is given in 
appendix A. From equation (A15) 


% = 


1 


2 

cos Y 


1 

TT * — ^ 


1 + ^ 


tan“^ 2H 



-b/2 


(15) 


where y is the dihedral angle in radians. Values of e for various H values 
are listed in table 2. 


End-Plate Wing 


An accurate but simplified e function for 
end plate wings is developed in appendix B. 

From equation (B14) 

2 -- n 

■ 0612 H(.8 - 
1 + .3688 


= 1 + 


"4 

TT 


1 + .732 H - 



H 


~r 

0 


b/2 


(16) 


Values of e for the end-plate wing computed from equation (16) are presented 
in table 2. 


90° Cruci form Wing and Circle Dumbbell Wing 

The 90° cruciform wing has H = 1, and as 
shown in appendix C, e = 2 for the unbanked wing. 

Also shown in appendix C is the variation of e 
with bank angle, and that when induced drag is taken as relative to wing length 
(and not lateral span) the e is unity at all bank angles. 

The e for the circle dumbbell wing can be developed from induced drag 
methods presented in reference 2. 
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GENERALIZED EQUATIONS FOR ESTIMATING e 


The objective is to develop an e function, with continuity, which becomes 
the proper e function of each of the basic configurations when that config- 
uration geometry is inserted. Because there are many basic configurations and 
the e function is continuous, the application to general configurations should 
be within an acceptable order of accuracy. The general e function is statis- 
tically developed so that each basic configuration is taken into account. 


e 



16 

2 


na H -)(1 - ^ 




(1 - .19 n-|Q H)(l + .4826 n^Q H^) 

(1 + 1.81 n-,0 H)(l + .255 ) 


(1 + f) + 


(17) 

(18) 


(6 - 4/2) (e, - 1) (e„ - 1) 


’af 




1 - /2 v/e 


- 1 - {/e~^^ 
V ' e 


2/e7^ 


-l2 


'af 


^ar 


naf=l/2 

1 = e^ - 1 


(19) 


where e^ and e^ are given in equations (15) and (16) and in table 2. 


f = 
P = 


1.4512 X I 


1 +3(1 - n 


p|.841 _ ,p, 1/2^1.0795 |-^ 

*^10 H 13/2 , 

r: FTaT2 (<^ 


+ .585|P|^( |P| - .37)^] (20) 

- a^) (21) 
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perimeter of full or of partial span symmetric ellipse boxwinq~ ] 

twice the wing span 


E. 


"i '’1 


HE , !L , , 

^ n- n- 


n 


1 


n 


1 


Also functionally 

„ (!L, , H_ (i.) ■ 

e'n-]^ n-| H 

Numerical values of Og/n-j are presented in table 3. 


An approximation for u is 


'1 


1 + 


(H_)1.11374 

"1 


1 


1 + ,H ^1.11374 


1.11374 


(^)'n 


1 + 


n 


n 


32 (1 + --)^ 
’^1 


For a semi-ellipse in a n-i-tiQ interval on the wing semi -span 

k 


ae = no + (n-, - Hq) 


E, 


Vi - 


ni - no 


ni - no 


,, _ cross-section area of boxwing 

M _ ^ 


HE . k 


ni - no 


H = ^ , where h^ is the maximum height of wing configuration, 


I (22) 


(23) 


(24) 


(25) 

(26) 
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( 27 ) 


n-|o = n-| - Hq, where n-| and ng are the spanwise stations of the 
outboard and inboard sides of the boxwing at the wing. 


n 


a 


,1 

nd(^) , averaged n of outboard face of wing. 

0 


(28) 


a = where p is perimeter of full or of partial span boxwing 
in which h is constant or decreases with n- 


(29) 


For bow-tie and dumbbell type wing configurations, the effective perimeter 
is determined from the shape with the maximum dumbbell height translated one- 
half dumbbell width towards the inboard side but limited to the inboard side 
of the dumbbell which is ng- Some examples are 


1 


1 

where- --is the translated shape needed to determine p and a for this type 

of wing. The value of of the dumbbell wing with symmetrical ellipses is 
given in equation (24). This and the dumbbell type wing value of a are 
used in equation (21) for evaluating P. 




n ^ = n , for full or for partial span boxwings in which h 

dT d 

is constant or decreases with n- 


\f 


1 + .05 


1 + .05 


^ 2 
"^me 


2 

^me 


for bow-tie, dumbbell type, and 
vertical fin wing configurations, 
where , 


(30) 


(31) 
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For fin wings 


, n is mean lateral span station 
for S-j = S 2 . For a vertical 
fin wing, = ^0 " ’^l " "^a‘ 


\f 


+ 2 (2 ig2 - Hq! 



2 

H 2 + (n-,2 - Oqi 


)‘ 


^'^11 ■ 


(1 


sin Y 2 
sin V-i^ 




.05 + 


(n-,1 - TiQi) hn 

‘~“~77r527 



+ 


1 .77 + .9H ^ 

9 1 + .67H' ^"^IS 



(nil + n ^2 “ ^’^Ol 


^ 

1 + .05 


2 

•^01 


+ 


.77 + .9H 
1 + .67H 



^2n-l " ^2n 


1.15 N H 


1 


(1 


^2n 

^2n-l 



(32) 


where the variables are defined through the following fin wing semispan 




TABLE 2. - IDEAL e OF V-WING AND END-PLATE WING 


H 

e 

V 

V-wing 
eq. (15) 

e 

e 

end-plate 

eq. (16) 

H 

e 

V 

V-wing 
eq. (15) 

end-plate 
eq. (16) 

0 

1 

1 

.8 

1 .32715 

2.38527 

.05 

1.00189 

1.09878 

1.0 

1.45265 

2.69328 

.10 

1.00750 

1.19523 

1.2 

1.58206 

2.99408 

.15 

1.01665 

1.28888 

1.4 

1.71362 

3.28952 

.20 

1.02909 

1.38192 

1.6 

1.84636 

3.58082 

.25 

1.04452 

1.47248 

1.8 

1 .97974 

3.86878 

.3 

1.06257 

1.56140 

2 

2.11345 

4.15394 

.4 

1.10524 

1.73481 

3 

2.78235 

5.5489 

.5 

1.15470 

1 .90313 

4 

3.44833 

6.9085 

.6 

1.20903 

2.06720 

CO 

-^2H/tt 

->4H/t7 


TABLE 3. - PERIMETER PARAMETER OF THE ELLIPSE BOXWING 


H 

'^l 

°e 

'^l 

H 

’^l 

o 

e 

"1 

H 

"1 

°e 

''I 

H 

''l 

°e 

"1 

H 

’^l 

'^e 

"1 

H 

^1 

°e 

'^l 

0 

1.000000 

.20 

1.050502 

.40 

1.150656 

.60 

1.276350 

.80 

1.418083 

1.0 

1.570796 

.01 

1.000275 

.21 

1.054653 

.41 

1.156442 

.61 

1.283110 

.81 

1.425490 

1.2 

1.73145 

.02 

1.000960 

.22 

1.058912 

.42 

1.162288 

.62 

1.289908 

.82 

1.432922 

1.4 

1 . 89808 

.03 

1.001977 

.23 

1.063276 

.43 

1 .168192 

.63 

1.296744 

.83 

1.440381 

1.5 

1.98318 

.04 

1.003286 

.24 

1.067740 

.44 

1.174155 

.64 

1 .303617 

.84 

1.447866 

1.6 

2.06931 

.05 

1.004856 

.25 

1.072303 

.45 

1.180724 

.65 

1.310526 

.85 

1.455376 

1.8 

2.24423 

.06 

1.006667 

.26 

1.076960 

.46 

1.186245 

.66 

1.317471 

.86 

1.462910 

2.0 

2.42211 

.07 

1 . 008700 

.27 

1.081 709 

.47 

1.192371 

.67 

1.324450 

.87 

1.470470 

2.5 

2.54309 

.08 

1.010940 

.28 

1.086547 

.48 

1 . 1 98549 

.68 

1.331464 

.88 

1 .478053 

3.0 

3.34118 

.09 

1.013375 

.29 

1.091470 

.49 

1.204778 

.69 

1.338512 

.89 

1.485660 

3.5 

3.81272 

.10 

1.015994 

.30 

1.096478 

.50 

1.211056 

.70 

1.345592 

.90 

1 .493290 

4 

4.28921 

.11 

1.018786 

.31 

1.101567 

.51 

1 .217383 

.71 

1.352705 

.91 

1.500943 

5 

5.25251 

.12 

1.021742 

.32 

1.106734 

.52 

1 .223758 

.72 

1.359850 

.92 

1.508619 

6 

6.22502 

.13 

1.024856 

.33 

1.111977 

.53 

1.230179 

.73 

1.367025 

.93 

1.516317 

8 

8.18624 

.14 

1.028119 

.34 

1.119741 

.54 

1.236658 

.74 

1.374232 

.94 

1.524037 

10 

10.15994 

.15 

1.031525 

.35 

1.122685 

.55 

1.243157 

.75 

1.381468 

.95 

1 .531 778 

20 

20.09713 

.16 

1.035068 

.36 

1.128146 

.56 

1 .249712 

.76 

1.388734 

.96 

1.539541 

30 

30.07148 

.17 

1.038742 

.37 

1.133675 

.57 

1 .256310 

.77 

1.396029 

.97 

1.547324 

40 

40.05720 

.18 

1.042543 

.38 

1 .139271 

.58 

1.262949 

.78 

1.403353 

.98 

1.555128 

50 

50.04800 

.19 

1 . 046464 

.39 

1.144932 

.59 

1.269630 

.79 

1.410704 

.99 

1 .562952 

00 

- H/n-, 
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APPLICATION TO PARTIAL SPAN BOXWINGS 


When a strut is designed so that it contributes to the lift of the vehicle 
and since the strut is attached to the wing, the configuration is the same as 
that of the boxwing, in particularly, partial span boxwings. One objective of 
this section is to demonstrate usage of the generalized equation (1 7) for 
estimating e for this type of wing. 

For the rectangle boxwing, = H, tiq = 0, n-| = n-|Q = = 1 , a = 1 + H , 

and equations (17) and (21) reduce to 


e = (1 + - H)(l + f) 

TT 

P = a - a 

e 


(33) 


Equation (33) is the same as equation (10) which is a solution for the 
rectangle boxwing. 


For the ellipse boxwing, H, = ttH/ 4, n, =irH/4, a = a„, then equation (17) 

da 6 

becomes e = 1 + H, which is equal to equation (6). 


For the diamond boxwing, = 1/2, = H/2 , Hq = 0, n-j = 1» n-]Q = 

and a = /I + H^, then equations (17) with (18) and (21) reduce to 


e 


fl + - H (1 - .44604 

L TT 


1 - .19H 1 + .4826H^ 

^ 1 + .255h5/2^J^ ^ 


(34) 


P = .35355 (/I + - Og) 


(35) 


where with the P of equation (35), f is determined from equation (20). These 
f values are negative since P is negative. Values of e computed from equation 
(34) compare within one-tenth of one percent accuracy with those of the dia- 
mond boxwing listed in table 1. 
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Reflexed Ellipse Boxwing 


For this wing, = Hgf = 1 - ir/4, 
Hg/H = 1 - tt/ 4, n-|g = 1 » a = then 
P = 0, f = 0, and equation (17) becomes 

e = , + (i- 1)H [l - (5 - f -f)sj 

= 1 + .27324H (1 - .881725 Gq) 



where Gq is given in equation (18) with ri-|Q = 1. Some typical values of e are 


H 

0 

.1 

.2 

.5 

1 

2 

e 

1 

1.00721 

1.02008 

1.07519 

1.19110 1 

1.46881 


Partial Span Ellipse Boxwing 



For this wing, = ^/4)n-| > H^/H = (Tr/4)n-] , 

= is a specified value, 0 s n-| ^1, 

Og, then P = 0, f = 0, and equation (17) becomes 



e 


1 + n-,H[l - (1 - n2)(l 


2 

16 ^1 



(37) 


When H = , the central wing is circular. For = 0.5, the equation and some 

values of e are 


e = 1 + .5H (1 - .63434 Gg) 
where Gq is given in equation (18) with n^g = 0.5. 


H 1 

0 

.1 

.2 

.5 

1 

1 

2 

e 

1 

1 

1.02116 

1.04710 

1.14372 

1.33842 

1 . 78399 
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Partial Span Semi-Ellipse Boxwing 


For this wing, = (TT/4)ni, 

(Tr/4)rii, n-jQ = n-| , then equations (17) and (21) 
become 



e = S 


1 + n 


1 H [l - (1 - ni)(l - 4 


P = 


4 *^1 


2^2 


1 + 3(1 - Jn^)(l - f n^) 


3/2 


(o - 




(38) 


(39) 


The semi-ellipse perimeter parameter is 







— 


1 + 

E, 

k - aJ 

TT^ip, 

"^1 

H ^ 

> 


m 

k-^ 


"^1 

, H 2 4- 



(40) 


and Og is given in equation (24). With n-j = 0.5, equations (38) through (40) 
reduce to 



(1 - .63434 Gq)J (1 + f) 

P = .033507 (a - a^) 


E, k = /I - (4H)Z, H <l 

O = T + T 

4 4 



4HE , k = A - (4H)-2, H > ^ 


(41) 
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Using equation (41) with equations (18) and (22), for n-| = 0.5, some typical 
values for e are 


H 

0 

.1 

.2 

.5 

1 1 

2 

e 

1 

1.02324 

1.05142 

1.15340 

1.35489 1 

1 .81216 


For the full -span semi-ellipse boxwing, i 
through (40) reduce to 

e = (1 + H)(l + f) 

P = .523103 (o - a ) 
e 

E, k = /I - (2H)^, H < 1 

a = i + Jr 

^ ^ 2HE, k = /I - (2H)-^, H > j 

Using equation (42), for n-| = 1 , some typical values of e are 


H 

0 

.1 

.2 

.5 

1 

2 

e 

1 

1.11669 

1.23959 

1.61191 

2.22938 

3.45462 

1 


a-j = 1 , then equations (38) 




(42) 


J 


Partial Span Rectangle Boxwing 


For this wing, = n-| , H^/H = o-] , 

a = n-| + H, n-iQ = n-| , then equations (17) and 
(21 ) become 



f 

-1 

+ f) 


1 


^1 

(43) 


( 44 ) 
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With n-i = 0>5, equations (43) and (44) reduce to 
e = [l + ^ (1 - .44604 Gq)] (1 + f) 

> 

P = .064 (H + .5 - Og) 


(45) 


where a is determined from equation (22) with n-| = 0.5. For n-| = 0.5, 
some typical values of e are 


H 

0 

.1 

.2 

.5 

1 

2 

e 

1 

1.05343 

1.10433 

1.26529 

1.55646 

2 . 1 8486 


Partial Span Diamond Boxwing 


For this wing, =(l/2)n-|. 

Hg/H =(l/2jn-|, a =/n2 + H2, n-jg = n-] > 
then equations (17) and (21) become 


2n-|H 




^ ^1 


n?)0 - I 


ni)G 


0 



P = 


1 2 
2 ^1 

1 + 3(1 - 1 ni)(l - 


-1 3/2 


2^2 


(/n^ + H2 



(47) 


With n-] =0.5 


(48) 

(49) 


H 

0 

.1 

.2 

.5 

1 

2 

e 

1 

1.00645 

1.01685 

1.06521 

1.17650 

1.44816 


e = [l + ^ (1 - .84251 Gq)] (1 + f) 

P = .011339 (/.25 + H2 - a^) 

For n-j = 0.5, some typical values of e are 



Partial Span Right Triangle Boxwing 


For this wing, =0/2)n^, ^ 

Hg/H = (1 H) , a = ( 1 / 2 ; m + ^1/4) + H2, -’l 1 

riiQ = n-i . then equations (17) and (21) 

become 



For n-j = 0.5, typical values of e are 


H 

0 

.1 

.2 

.5 

1 

2 

e 

1 

1.00699 j 

1.01853 

1.07080 

1.18579 

1.46190 


For = 1 > typical values of e are 


H 

0 

.1 

.2 

.5 

1 

2 

e 

1 

1.03178 

1 .07181 

1.23909 j 

1.60153 

2.37390 


Ellipse Dumbbell Boxwing 


In this wing 

^a " ^af " I (^1 ^"^1 ■ ^0^ 

Ha 

. H“ " 4 ('^1 ■ '^0^ ’ "^lO " "^1 ■ ’^0 

and a=ag, soP=f=0 
where is from equation (24). 
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With these values and using equations (17) through (19), e can be 
determined for arbitrary ellipse and position in the wing semispan. For the 
special case of circular boxwings with n-| = 1 » the wing appears as shown below. 


O ^ Q-—0 

% 



For the circle dumbbell boxwing with n-j = 1 , the height parameter is 
H = (1 - t]q)/Z, and the wing terms simplify to 

(1 + tiq) + t (1 - ng) 


9 _ T7 


H 4 " ’^0^ 


(53) 


nio - 1 - Iq 


riio H 2 '^0^ 


With the values of equation (53) and equations (17) through (19) typical 
values of e are determined at hq = -8, .7, .5, and 0, which correspond to H 
values of .1, .15, .25, and .50, respectively. 


H 

0 

.10 

.15 

i 

.25 

.50 

e 

1 

1 . 1 7494 

1.24393 

1.39304 

1.64077 

ref. 2, e 

1 

1.175 

1.254 

1.393 

1.630 


These circular dumbbell boxwing values compare well with the theory results 
from the application of the methods of reference 2. 
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Rectangle Dumbbell Boxwing 


In this wing, 

"a = '^af = 

H 

'^10 " '^1 ■ '^0 = yr 

a = !>! + H 




^0 

n 






1“ 

1 





(54) 


With the terms of equation (54), e is determined by using equations (17) through 
(21) and (24). For the special case of a square dumbbell boxwing with = 1 , 
then H = (1 - riQ)/2 and the equations simplify to 


e = + f [l - (1 - ^H)(l - 2 H)Gq] 

[l - (2H)^'^] (e^ - 1) 

G = (1 - -38 H^)(l + 1,9304 H^) 

° (1 + 3.62 H^)(l + 1.4425 H®) 

a = 1 + H, a = 1 + .4221 T2H 
e 

p ^ 8(.577888)H^ 

[l +3(1-4 H^)(l - 8 

where H = (1 - nQ)/2. With equation (55), typical values of e of the square 
dumbbell boxwing are determined at ng = .8, .7, .5, and 0, which correspond to 
H values of .1, .15, .25, and .50, respectively. 


H 

0 

.10 

.15 

.25 

.50 

e 

1 

1.20451 

1.29238 

1.50320 

2.00026 


(1 + f) + 
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Bow-Tie Boxwing 


For this wing. 


H. 

" ’^af " ^ ’ FT 


1(1 


no)(l -fr) 


'10 


1 - Hq* cr - ^ [l + Hq ■*■ H + 


1q + VCh - Hq)2 + (1 + ■ 



(56) 


With the terms of equation (56), e can be 
through (21) and (24). For the special 
case of Hq = 0, and Hq/H = (tt^ - 8)/8 = 
.2337, the equations for determining e 
s i mp 1 i fy to 


determined by using equations (17) 



e 


P 


(1 + J H)(l + f) + .560149 (e^ - 1) 


E, 

HE, 


,44219 




+ .61685 H + 


^/r 


+ .58722 HZ 


(57) 

k = 7l - h2 1 

k = i/l - H-2 J 


Typical values of e for this bow-tie boxwing computed from equation (57) are 


H 

0 

.1 ^ 

.2 

.5 

1 

2 

e 

1 

1 . 24446 

1.45077 

2.03665 

2.94646 

4.69481 


APPLICATION TO FIN-WINGS 


Fin-wing configurations include end-plate wings, winglets, V-wings, cross- 
wings, and fan wings. Since in this class of wings H = 0, equation (17) 
simplifies to 

e = 1 + G.| (58) 
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where G-| as a function of is given in equation (19), and in equation 
(32). The prediction of e for fin-wings can be demonstrated by application of 
the method for estimation of 


End-Plate Wings and V-Wings 


For the end-plate wing 

'^al " '^a2 " '^ll " '^01 " '^13 " ^ ’ 

Y2 = -Yi, 

then equation (32) reduces to 

H2 

1 + .05 

^1 

'^af " 1.05 , at H 2 = H.| , - 1 



At n 


af 


1 , G.| = e^ - 1 , so by equation (58), e = e^, (eq. 16). 


For the V-wing 


■^al "^a2 -5’ n-|^ 1 , ngi n-jg = 

then equation (32) reduces to n ^ = 1/2. 
so by equation (58), e = e^, (eq. 15). 

Cruci form or 

For the cross-wing 

’^al " "^a2 " '^11 " ^^12 " ^ ’ 

■^Ol " '^13 " 

then equation (32) reduces to 


0 » Y2 Y"! » 

At this value of ^1 ~ ^v ~ ^ ’ 
Cross-Wing 
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^ [l + .257267 


(60) 


/lZ3 

Vl + 


0 


sin Yq 1 


When Y 2 ~ Y-] > then = 1/2, the configuration becomes the V-wing, and as 

before, e = e^. When Y 2 ^ -Y^l and H 2 = H-j , the configuration becomes a 
symmetrical cross-wing, and equation (60) leads to = .757267. With this 
value of and equation (19), equation (58) gives e for the symmetrical 
cross-wing as 


e 


1.84922 (e^ - 1) 
(.397872 +V^-~T 


(61) 


where e^ and e^ are given in equations (15) and (16), and in table 2. Some 
typical values of e are 


iL 

0 

.1 

.2 

.5 

1 

2 

e 

1 

1.03997 

1.11847 

1.43415 

2.00000 

3.09219 


where the H = 1 configuration is a 90 degree symmetrical cross-wing which as 
shown in appendix D has a value of e = 2. 


For this wing. 


Winglet Wings 
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then with equation (32) 


n =[l + 1 n2 (1 _ !!!L!2) 

V 2 '^ 2 y„ 2 ,(, . slny,) 


(1 - loi^ ■ *^01 I 


2.1527 


' ^ "01 

2(1 + .05 Hq,) 


For a single winglet, y 2 " > then 

'' ^1 

-f ? 

2(1 + .05 Hot) 


At H = 1, equation (58), with (19) an<i the e^ and e^ values at H = 1 , results in 


e = 1 + 


2.14961 Hgf 
(1 + ,12672Vf;^^)‘ 


Typical values of e for various locations of a single winglet for H = 1, are 


'^01 

0 

1 .25 

.50 

.75 1 

1.00 

1.05 

■^af 

. 50000 

.62305 

.74074 

.85106 

.95238 

j 

.97145 

e 

1.45265 

1.68961 

1.95892 

2.24811 

2.54421 

2.60315 


where the Hq^ = 0 value is the same as for the V-wing. 


For vertically symmetric winglets. 


Y 2 ^ 1 * ^12 ’ ^2 

equation (62) becomes 


1 , H 2 = H.| , then 



1.05 + 


(1 - nQ^) n - ngi 


2.1527 


I 

■ 2M + 


1 + 


2(1 + .05 non) 
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Then using equation (64), typical values of e for various r^-j values of a wing 
with vertically symmetric winglets, for H = 1, are 


'\)1 

0 

.25 

.50 

.75 

1.00 

1.05 

^af 

.757267 

.84219 

.89943 

.94303 

1.00 

1.01498 

e 

2.00000 

2.22356 

r 2 . 38588 

i 

2.51568 

2.69328 

2.74146 


where the Ogi = 0 value is the same as for the symmetric cross-wing of H = 1, 
and the ngi = 1 value is the same as the end-plate wing of H = 1 . 


For the case of y 2 = "^IZ ^ 
then ^2 = 0, and equation (62) becomes 


\f 


i / ~ *^01^ 

'Vh 2 Td - ^)2 


.05 + 



1 + 


'01 


2(1 + .05 ngi) 


( 66 ) 


Then using equation (64), typical values of e for various ngi values of a wing 
with upper surface winglets, for H = H-j = 1, are 


^01 

1 

0 

.25 

.50 

.75 1 

1.00 

’ 

fi) 

1 j 

.59096 

. 68880 

.77623 

. 86222 

.95238 

e j 

1.62335 

1.83499 

2.04810 

2.27930 

2.54421 
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Fan Wings 


To analyze a 2N - 1 fan wing where 

'^,2n-l " *^,2n " ’iD.n-l " "bl ’ 
"^2n" ~^2n-r ^2n-l " ^2n ’ ’^a,2n " 
^a,2n-1 -^1/2X1 

then from equation (32) 



"af 



(1 - Hq^) 

2.1527 


3/2 


+ 


2 .77 + .9H 
9 1 + .67H 


0 



1 + 

2 + .1 



+ 


.77 + .9H ^ 

1 + .67H 1.15 N 



(1 


.36 ^ 

2.1527 -J Z.3 


’2n-l 



(67) 


For equally spaced fins along H 


^2n-l _ n - 2 
H.| N - 1 



'2n-l 



= 1 (N - 2) 


(68) 


For a three-fin fan wing, N = 2, then from the first part of equation (67), 
the values for various ng-j span stations and typical e values at H = 1 
from the equation (64) are 


"01 

0 

.25 

.50 

.75 

1 

1 

*^af ! 

.86838 

.92008 

.94059 

.95486 

1 

e 

2.29667 

2.44668 

1 2.50828 

2.55182 

j 2.69328 
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For other values of H , e can be determined from equation (58) with G-j given in 
equation (19), where are those listed above. 

For a five-fin fan wing, N = 3, and with equations (67) and (68), the n^f 
values and typical e values at H = 1 are 


"^01 

0 

.25 

.50 

.75 

1 

^af 

.94295 

.98166 

.98352 

.97640 

1 

e 

2.51544 

2.63513 

2.64098 

2.61862 

2.69328 


For an in finite- fin fan wing, N = <», then with H = 1 


^1 

0 

.25 

1 

.50 

.75 

1 

"af 

1 .0921 

1.1048 

1 .0694 

1.0195 

1 

e 

2.9992 

3.0433 

2.9216 

2.7561 

' 

2.6933 


The = 0 value of e is the same as the bow-tie wing with ng = Hq = 0, 
and H = 1. All the above r^-] = 1 values of e correlate with the e of the 
symmetric end-plate wing with H = 1. 

For the infinite-fin wing with dqi = 0, equation (67) simplifies to 

.75727 + .3348 (69) 

Then with various H values, typical values of e from equations (58) and (19), 
are 


H 

0 

.1 

.2 

.5 

1 

2 

\f 

1.01507 

1.02712 

1 .03775 

1.06324 

1.09207 

1.12498 

e 

1 

1.2525 j 

1 .4575 

2 . 06.77 

2.9992 1 

4.7419 
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These values of e for the infinite-fin wing with = 0, are close to those 
listed for the bow-tie boxwing where Hq/H = .2337, hq ^ is apparent 

cross-wing have the same e. 


the Hq = = 0 bow-tie and infinite- fin 

An example usage of the fin wing 
equations for other shapes is the case 
where 

*^11 ~ *^1 2 ~ ^ ’ ^2 ~ " Y ] 5 ^2 - H 
^01 " ^13 " ^15 " ^16 " '^02 

’^6 "" “'^5 " ^5 ^6 "" ^1 ’ 

\5 " ’^ae " 

then equation (32) reduces to 
n - 75727 + 3344 



(70) 


Equation (70) is the same as equation (69), hence e is equal to that listed 
under equation (69) . 


APPLICATION TO MULTI-WINGS 


Multi-Wings with Equal Spans 

A vertical stack of wings with equal spans is pictured as 


-b/2 b/2 

I T — 



r ‘ 

> 

V • 1 

wing 1 
' wing 2 

i 

*^3 



' ""wing 3 

3 

^ ^N . 

'''^wing 4 

f 


wing N 
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then = H, n, = 1, and equations (17) and (21) become 

a a 


e = (1 + ” H)(l + f) 


(71) 


P = 


(72) 


Comparing equations (71) and (10) shows, for small H so that P and f are small, 
that the e for multi-wings of equal spans approaches that of the rectangle 
boxwing. This is in accord with multi-plane thoery. The perimeter p for 
multi-wings requires a redefinition here since there are no vertical panels. 

In effect the h distances in p are subdued so that for equal span wings 


= 1 + 


1+3.06 


H o 
(Jl)2 


(73) 


where denotes a of a multi -wing aircraft and N is the number of wings and 
N - 1 is the number of spaces between wings. For equally vertically spaced 


wings , h^ = h^ 


hj^j, then 


^(^)2. (, _ i)(^)2._^_^ 

^ ^ (N - 1)^ 


r ^2 12 _ 

Ih7(n - i)J 


N - 1 


thus for wings of equal spans which are vertically equally spaced, the effective 
perimeter parameter is 




1 + 


H 

3.06 .1/6 
N ^ 


(74) 


The e is determined by inserting of equation (73) or (74) into equation (72) 

and using equations (20) and (71). Equation (71) is accurate in the range 
0 £H < 1, but somewhat higher as N becomes large. As mathematically proven in 
appendix D, in the limit for H very large, the e of the multi -wing and the 
lifts of wings of equal spans are simply 


29 


(75) 


e = N, and -r^ = 1 

which is independent of distribution. As the number N of wings becomes 
large, then as shown in equation (74), a= 1 + H, which is also the value for 
the rectangle boxwing. Thus the e of the multi-wing with an infinite number 
of equal span wings is identical to the rectangle boxwing of the same H. 

Multi -Wings with Arbitrary Spans 
A vertical stack of wings with .various spans is pictured as 



h2 

r , 

r 


^3: 



^ : 

r 


r 5 

i ' 

^4 

r 


E .1 ^ 

(76) 


The effective integrated value for n 

a 

is the sum of the rectangle areas formed 
by each of the multiwings, then 



where H^/H = (h^/b)/(hyb) , and b is the span of the widest wing. The 

|r 


effective perimeter of the multiwing 
is approximated by 


perimeter 


H2 + 


+ H2 + 






n+l 


(77) 
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The multiwing perimeter parameter is 




ytcr - r\2)^ + (1 - - (1 - no)B 


Ti2> 

i + B 


ri2> 


(78) 


where 

(2 - nn)^(^)^ (79) 

As N approaches infinity the multi-wings outline a boxwing and equations (76) 
and (77) become the values of the boxwings. Then since at N = «, B = 0, the 
prediction of e for the multi-wing of infinite wings is identical to that of 
the boxwing formed by the outline of the infinite wings. 

The n-jQ parameter is ri 2 which represents the span extent of the wing 
with second largest span, then 


B = 3.06 


N - 1 ul/6 


nio H2 


(80) 


The elliptic parameter is for a partial span to the ri 2 or n-|Q span station. 
The values for are determined from using a^/n-iQ instead of Og/n-j in 
equation (22) or in table 2. Thus with q and H /H from equation (76), and 

a d 

from equation (78), then P is evaluated from equation (21) which becomes 




'^lO'^a 


3(1 - n^on,)(l 




13/2 


■ "^10 


'^lO 


(81) 


With P, f is found from equation (20), then e computed from equations (17) and 
(18). 


Accuracy is ensured to H values such that e/(l + f) is the same as the 
H = «> value of e. The H = value is derived in appendix D where from equation 
(D9), for maximum e at H = “ 


e 




(82) 
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For finite height parameter H, e and are less than these values. 
The biplane solution has N = 2, then 


^2 ~ ’^a " ^ a ^^ ~ ’^10 ~ ~ ^^2 

equations (77), (78), and (79) become 



a = h 2 + H , B = 1.53(2 - 


" ^2 ^ 


.pTo 


ri2)^ - (1 - 0^)5 

nn 


(83) 


With these values P is determined from equation (81) and e from equations (17), 
(18), and (20). 


Diamond Shaped Multi -Wings 

A diamond shaped multi -wing can be 
represented by letting wing gaps and 
wing spans be >— 



_ (1/2XN + 1) + 1 - n I 

^n a/2j(N + 1) J 


» rj r 

_ ’'3 ^ 

’^.5(N+1) 


(84) 


where N is an odd integer. Then from equation (76) 


H 


(1/^N+l) 


n=2 


H 


% r " ::2 


4 0/iNtl) „ , 3 


N - 1 


n=i 


- n) = 



[<' 


+ 1 


- 1 )(- 


+ 3^ 


1 /N + Is/N + 1 

2 ' 2 2 



(85) 
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(86) 


or for all odd integer N, = H^/H = 1/2. The ri^Q value is 

_ N - 1 
"^10 "^2 N + 1 

The perimeter is 

° ~ ’^N + 1 ^ ^^2 ^ •\/^n2 “ nj^j 4- i) + (H - 

~T~ 


2 2H 

N + 1 N - 1 


+ 



N - 3^2 
N + r 


+ ( 


N - 3x2 
N - r 




(87) 


The multi-wing perimeter parameter is determined from equation (78) in which 
^2 given by equation (86). B from equation (79) summed for this diamond 
shape is 


_ 6.12 1 + f7/$i 

N 1 + N 


( 88 ) 


With aj^ evaluated, P is computed from equation (81) in which = 1/2, n-jg is 
from equation (86), and is with the ellipse extending to n-io* 

When N becomes infinite, these parameters become n = H /H = 1/2 (for any 

S d 

N) , n-|Q = 02 “ i > 8 = 0 , = o = /I + H2 . These values are the same as those 

listed in equations (34) and (35) for the diamond boxwing, hence e is the same 
as those listed in table 1 for the diamond. 


Numerical Values of e for Multi -Wings 

For equal span wings and vertically equally spaced wings, for the biplane 
N = 2, and for the triplane N = 3, equation (74) becomes 




H 


1 + 1.53 H 


yg (biplane) 


^ 1 + 


H 


1 + 1 .02 H 


jjQ (tri plane) 


(89) 
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With equation (89), P is evaluated from equation (72) with of the ellipse 
extending to n-j = 1. With P, f is evaluated from equation (20), then e from 
equation (71). These e's are listed in table 4 for various H values with the 
H = “ values from equation (75). 

For a partial span wing on the biplane, N = 2, = H, = H^/H = = 

r\^, and is given in equation (83), where denotes the spanwise extent of 
the lower wing. Equation (81) becomes 


P = 


LI + 3(1 - 


no) (1 


24)' 


■] 


3/2 


(of/i - 


ae) 


(90) 


where is with the ellipse extending to the n^Q = span station. With P, 
f is determined from equation (20). Equation (17) becomes 


e 


"tv [' - (' -T''2)C - 

TT 



(91) 


where Gq is given in equation (18) in which n^Q = ^ 2 - The e for the biplane 
with ri 2 = 0-6 computed from equation (91) is listed in table 4 for various 
H values. The H = “ value of e is from equation (82). 


TABLE 4. - IDEAL e OF VARIOUS MULTI-WINGS 


H 

equal span 
bi plane 
eq . ( 89 ) 

equal span 
tri plane 
eq . ( 89 ) 

ri 2 = 0.6 span 
^ biplane 
eq. (91) 

0 

1 

1 

1 

.1 

1.2021 

1.2176 

1.0456 

.2 

1 . 3529 

1 . 3834 

1.0876 

.3 

1.4810 

1.5315 

1.1301 

.4 

1.5842 

1.6654 

1 . 1 769 

.5 

1.6407 

1.7839 

1.2154 

oo 

2 

3 

1.36 
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CONCLUSIONS 


It has been demonstrated that the application of the generalized equations 
(17) through (32) and equations (76) through (82) provides a direct method for 
estimates of ideal e for arbitrary lateral-verti cal wing forms. Example pre- 
dictions of e are made for a wide range of wing configurations, from multi-wings 
to fan and dumbbell wings. Several check points were verified with other 
methods including the cruciform at H = 1 (eq. 61), the circle dumbbell (eq. 53), 
and biplane and triplane (table 4 with ref. 1). Also, multi-wing solutions 
converge to those boxwings as the number of wings increase, and fan wing solu- 
tions converge to those of bow-tie wings as the number of fins increase. For 
the V-wing, an exact closed function was developed. An asymtotic unique maxi- 
mization solution for multi-wings was found which resulted in exact expressions 
for ideal e and ideal lift of each wing. 

Hampton Technical Center, Kentron International, Inc. 

Hampton, Virginia 23666 
June 20, 1980 
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APPENDIX A 
V-WING e FACTOR 



The dimensionless added mass is given by 


M 1 - 


2 2 
PIT a a a 


where a is determined from the definite integral 
1 

fdg 


(1 - c)'-“(c + e)“ 


where 


_ 1 - a n 1 

C - — , 0 ^ a ^ ^ 


(A2) 


(A3) 


(A4) 


This integral was evaluated numerically in reference 3 for various values of y, 
then dimensionless added mass was computed and tabulated. However, as will be 
shown here, the integration in equation (A3) can be obtained by expanding the 
second parenthesis term into a binomial series. Thus 


(c + d' 




(a + 1) ± . . . 

2!c 


(ct+l)(oi+2) . . . (a+n-l)^*^^ + . . . 

n !c^ 


(A5) 


36 



The definite integrals of the following type can be evaluated: 


f r' 

d? 

n=0 

1 


n=l 

2! 


J (1 - 

^)l-a 


a( a + 

T7 

a( a + 

1 ) { a + 

■27 j 

0 









n=2 



3! 


{ 




a( a + 

T)T 

a + 2) ( a + 

■ 

( 



n=n 



(n + 1)! 

1 





a( a + 

Trr 

^T+ 2) . . 

. ( a + 

J 


Inserting equation (A5) into equation (A3) and applying equation (A6) 
® ^ [“«(a + 1) ■ l)(a + 2)c + 2)(„ + 3 )c2 


(-l)'^(n + 1) 

(a + n)(a + n + l)c*^ 


With c given in equation (A4) , equation (A7) can be written as 


. = ( 1 - V (-l)^(n + 1) , a ^n+^ 

a ( n + a) ( n + 1 + a) 1 - ct 


For the range of 0 ^ a ^ 1/2, this sunmation is identically unity, 
exact solution for a is simply 


a = (-^^) 

a 


With equation (A2) the dimensionless added mass is 


M 


pira 


2 M - a'' 


(A6) 


leads to 


(A7) 


(A8) 

Thus the 
(A9) 


(AlO) 
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From reference 3, the induced drag is given by 



u. 2 

^ 2 MV 

which equated to equation (1) results in 


From the V-wing sketch the geometrical a and y are 


(All) 


(A12) 


^ = 1 
b ' cosy 


(1+4 , Y = tan"^ 2H 


(A13) 


With equation (Al), the dimensionless added mass and factor e, in terms of the 
Y angle in radians, in closed functions, are 


M 

2 " 
pira 


1 


1 



cos Y 


2y] 


TT 



ly 

TT 


TT J 

1 - 
1 + 



TI 




(A14) 


(A15) 


In terms of the height parameter, H = h^/b 


pira 


1 - - tan"^ 2H 

TT 


1 + - tan"^ 2H 


2 -1 
- tan 2H 

TT 


e = (1 + 4 H ) 


2 -1 

1 - - tan 2H 

TT 


1 + - tan"^ 2H 

TT 


2 -1 
- tan 2H 

TT 


(A16) 


(A17) 


Equations (A15) and (A17) give e in terms of elementary functions of either the 
dihedral angle or the height parameters of the V-wing. 
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APPENDIX B 


END-PLATE WING e FACTOR 


From page 211 of reference 1, 
Y = 6k sn n, -| = 


then 


sn 



/ K - E v1/2 _ 


sin <#> 


(B1) 


(B2) 


From the definition of snn related to sin<j>, n is 


= F( (})\a) = 


(1 - 


.2 .2 , 1/2 

sin oSin ({)) 


(B3) 


that is, the incomplete elliptic integral of the first kind. In addition, 
from reference 1 , 



where k = sina, k' = (1 - = cosa. With equations (Bl), (B2), (B3) , the 

h^ of equation (B4) becomes 


H = ^ [E(<J>\a) - I F(<f\a)] 


(B5) 


(p = sin 


’ 1 
sina 



(B6) 


where F(<p\a) and E(<^\a) are incomplete elliptic integrals of the first and 
second kinds with amplitude and modular angle a, and K and E are complete 
elliptic integrals of first and second kind with modular angle a. The term 
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within the brackets of equation (B5) is the Jacobian Zeta Function. Then 

H = |KZ(<f>\a) (B7) 

where the function KZ( (j>\a) is tabulated in reference 4 for various values of 
^ and a, and has a maximum value at the <l> given in equation (B6) . From 
equation (B4) 

4 |,2 /2E 2 V /Do^ 

e = K ("j^ - cos a) (B8) 

7T 

As can be seen from equations (B5) through (B8) , the factor e is not a 
direct function of H but is related to H only through the parameter a. That 
is, a value for a is specified, then ((> is determined from equation (B6) , then 
H from either equation (B5) or (B7), and e from equation (B8). Some typical 
results for a range of a values are shown as follows: 


a 

deg 

H 

e 

r 

a 

deg 

H 

e 

r 

30 

.07194 

1.14138 

.6996 

80 

.89782 

2.53705 

.7312 

45 

.1 7351 

1.33322 

.7082 

82 

1.0165 

2.71842 

.7296 

60 

. 34824 

1.64573 

.7200 

84 

1.1718 

2.95205 

.7263 

68 

.49550 

1.89566 

.7268 

85 

1.2713 

3.10002 

.7233 

75 

.68710 

2.20712 

.7318 

88 

1.7821 

3.84313 

.7048 


This list includes the parameter r which is evaluated from 

= 1 + 2H - e 

(e - 1 - ^ H)H 


(B9) 


Equations (B5) through (B8) are quite complicated functions and do not 
permit e predictions for a specified value of H. It has been found that a 
relatively simple function exists between e and H. In the limit as a ->■ 0 
equations (B7) and (B8) become 

e -y 1 + ^ 
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2 

where k = sina -^0. Eliminating k results in 
H->0 

e = 1 + 2H (BIO) 

For a ^ ir/2, equations (B7) and (B8) become 


H 



e 



4_ 

k' 


where k' 


(1 


k2)V2 = 


cos a. 


Eliminating the log term, e becomes 



7T 


(Bll) 


A combination of equations (BIO) and (Bll) indicates a function of the follow- 
ing type 


e = 1 + 




rH 


) H 


(B12) 


where r is a quasi -constant. The solution of equation (B12) for r is equation 
(B9). As can be seen from the listed H and r values, r is nearly constant 
with H, but remains a weak function of H. An approximation for r which dupli- 
cates those in the list is given by 


■0612 (.8 - H)^ 
1 + .3688 


(B13) 


Then equation (B12) gives e with a four to five digit accuracy for the wing with 
end plates as 


e 


+ 



1 + .732 H 



.0612 H(.8 - H)^ 
1 + .3688 


(B14) 
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APPENDIX C 


90° CRUCIFORM WING 


As discussed in reference 5, page 39, in a 90° cruci form wing there is no 


induced velocity at one vorticity sheet 
due to the other. In this case each 
wing acts independently of the other, 
thus for minimum induced drag the 
spanwise loading distribution is ellip- 
tical along b . With L . and L „ 
representing the lift normal to wings 
1 and 2, with equation (1) 




2 2 2 

since L^.| + 1^2 " *- • Thus for this case e = 1. In reference 5 when the 
C| , Cn 5 and A are reduced to dimensional quantities with b, ,, the solution is 
also e - 1 . L jL 

For a 90° cruciform wing banked L . / , \ 

^n2 

at angle <j>, the loading is elliptical 

along wings 1 and 2, and, similarly / 

to equation (Cl), the induced drag is wing 1/ wing 2 

0.=—^+^^=-^ (C2) 

irq b Tiq b irq b 

^ w ^ w ^ w 


wing 2 


2 2 2 

since again L^.| + L ^2 ^ L , and hence e = 1. Thus when the induced drag of a 

90° cruciform is in terms of the wing width, b , the efficiency factor is 

w 

unity at any bank angle. 

Generally the induced drag of a lifting configuration is in terms of the 
lateral span of the configuration, that is, b. Then from equation (1), the 
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I 


induced drag of the 90° cruciform wing is 


‘^i 2 — 

frq e 

Equating equations (C3) and (C2) results in 

bo 1 

e = ^ 1 _ 

cos2(f-U|) 

where (j> is the angle of bank of the 90° cruciform wing. When (j) = 0, 
±ir/4, then e = 2, 1.176, and 1, respectively. In summary for the 90 
cruciform wing 

e = 1, for all bank angles when is defined by equation (Cl), 

e = equation (C4), for all bank angles when D. is defined by 
equation (C3). ^ 


(C3) 

(C4) 

±ir/8. 
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APPENDIX D 


ASYMPTOTIC SOLUTION FOR OPTIMUM e AND LIFTS 
FOR MULTI-WINGS WITH LARGE GAPS 


If the gap height h^ of a multi -wing is large enough such that the 

influence of one wing on another is negligible, then the ideal loading for 

2 2 

each wing will be elliptic and induced drag is given by D^. ^ = L^/Tiqb . The 
sum of these induced drags equals that of the multi-wing having the sum of all 
the lifts. Thus for a multi -wing with N planar wings 




T:q b.j e n=l -nq b^ n=l 


(Dl) 


where b.| is the wing with maximum span. In dimensionless coordinates 

L b 

0 _ n _ n 

n L-j n b.| 


(D2) 


then equation (Dl) becomes 

(1 + i ^n) ^2 

e n=2 


solving for e 


(1 + i 0‘ 

^2 


e = 


1 + 


n=2 ri„ 


(D3) 


Equation (D3) is valid for any lift variation among the wings, however, there 

is some lift variation for which e is maximum. For equal loading e can be 

-2 - 2-22 

less than unity when 1 + ri 2 + • ’^m > • The maximization tech- 

nique of reference 5 can be applied to equation (D3). A partial derivative 
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of e with respect to each of the loading ratios gives 




Q ) jn 
2 


The conditions for maximum e are that these derivatives be zero for each of 
the loading ratios. Then 


N 

^2, 


- (1 


^2 


£ ) E 0 
r\' 2 

'^m 


(D4) 


These are non-linear simultaneous equations in with n = 2, 3, . . . N 
unknowns, and m = 2, 3, . . . N equations. For the biplane, N = 2, then 
equation (D4) with (D3) leads to 


1 (1 + 


2 

U2 






e = 1 + n 


(D5) 


For the triplane, N = 3, then the solution of two equations of (D4) leads to 


£ 


2 


1 


1 + £. 


2 

"2 


2 

02 





1 + 02 - 
2 “ 
^3 



(1 + 02) 



(D6) 


which involves the solution of a cubic equation. As can be seen, the solution 
of equation (D4) becomes complicated when N s 3. What is needed is to have 
£^ as some function of so that equation (D4) is satisfied. Such a function 
has been discovered by the hint given in equation (D5). Let 
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(D7) 


then equation (D4) becomes 


N 


2 ^ 2 . 


E 0 


(D8) 


n=2 


'm 


which completely satisfies the maximum conditions, and hence equation (D7) 
gives the individual wing lifts for maximum e. With equation (D7) inserted 
into equation (D3) the maximum or ideal e is 


i 


n 



e 


e 


max 



(D9) 


For equal span wings, = 1 , so 


i = 1 , e = N 
n 


For biplanes 


^2 ™ ^ 


= 1 + n' 


For triplanes 


^2 ^ 2 * 


2 2 2 
n^j e - 1 + rt 2 + ^3 


and so forth. 


(DIO) 


(Dll) 


(D12) 
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